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METHOD AND SYSTEM FOR DIGITAL PRINTING BY DENSITY 

SCHEDULING 

BACKGROUND OF THE INVENTION 

Field of the Invention 

The present invention generally relates to method and system for digital printing, and 
more particularly to a method and system for digital printing, mostly addressed at printing 
images, which exploit better the two-dimensional aspect of images than error diffusion methods, 
and take much more account of the image than dithering methods. 

Description of the Related Art 

Most printers today can print in only a limited number of colors. Digital half-toning is a 
technique for printing a picture (or more generally displaying it on some two-dimensional 
medium) using small dots with a limited number of colors such that it appears to consist of many 
colors when viewed from a proper distance. For example, a picture of black and white dots can 
appear to display gray colors when viewed from some distance. 

The fastest and most commonly used methods for digital half-toning are dithering 
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algorithms Which use threshold arrays (i.e., also called dither matrices or dither masks). 

The principle of this conventional method, as illustrated by Figure 1, is well known. The 
method allows to associate a matrix N of discrete values at 13 (typically "0" or "1", where "1" 
means a pixel is printed, and a "0" means nothing is printed) to an image I at 1 1 using a dithering 
mask 12 (i.e., a usually smaller matrix of threshold values). Various masks can be devised, 
according to the needs of precise applications, and several methods to devise masks with good 
performance have been disclosed (e.g., see U.S. Patent No. 5,1 1 1,310 to Parker et al., U.S. Patent 
No. 5,917,951 to Thompson et al, and U.S. Patent No. 6,025,930 to Thompson et al, each 
incorporated herein by reference). 

However, instead of a dithering mask, one can also use other half-toning algorithms such 
as an error diffusion algorithm. 

Before defining this class of techniques, it is noted that when using dithering, one could 
treat each pixel at any stage of the process, or treat a lot of the pixels in parallel, as the printing 
decision at each pixel (i.e., to print or not to print at that place) does not depend of the decision at 
any other point. In practice, one treats the pixels one-by-one, line-by-line, but this is just a matter 
of practicality, not a need of the method. 

To the contrary, error diffusion is rasterization-dependent. That is, one can choose an 

order in which the pixels will be dealt with, e.g., line-by-line going down the page, and from left 

to right on each line. Once this order is chosen, before making a printing decision at pixel (i, j) 

(i.e., choosing which color(s) will be printed at that pixel, if any), one modifies the input at (i, j) 

by weighted errors made previously in the neighborhood, and computes a new error as the 

difference between the modified input and the printing decision. 
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Dithering iriasks as well as error diffusion can be easily adapted to use more than two 
possible outputs per pixel (i.e., several levels of gray instead of just black and white). In such a 
case, one then speaks of multi-tone printers (with a digital printer then referring to black and 
white pixels only). Also, all these techniques, can be easily adapted to color printing. That is, 
both adaptation to multi-tone and to color (i.e., digital or multi-tone color) are well known by one 
versed in the art of digital printing. 

A well-known difficulty associated with dithering for a color image is linked to the fact 
that shifts on the paper, or the printing guns, may occur between different color components, 
which is known to create undesirable moire patterns. This problem is usually corrected by using 
different masks in the various color planes. 

Most of these techniques are reviewed in the book Digital Halftoning, MIT Press, 
Cambridge, MA (1987) by R. Ulichney, which is a general reference for digital half-toning in 
black and white. Regarding color, "Optical Color technology for electronic imaging devices" by 
Kang, Henry R. (Bellingham, Wash., USA: SPIE Optical, 1997) is a general reference for color 
digital printing. 

Other, more global, and very much more time-consuming, methods have been proposed 
to perform digital half-toning, such as the use of neural networks or global optimization. 

Thus, it is easy to understand why dithering works well for a uniform grey input. To the 

contrary, one can easily exhibit examples of images whose rendering by digital printing would be 

arbitrarily bad if using a dithering mask. This is because dithering only considers one pixel at a 

time. The reasonable results one obtains in practice for a natural image, such as pictures of 

natural scenes, are linked to the fact that, in general, such pictures are quite smooth. For 
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purposes of the present invention, "smooth" refers to a generally slow variation of the output 
from a pixel to the neighboring ones. Some poor cases are nevertheless commonly encountered. 

Error diffusion avoids this problem to some extent, but suffers from the fact that it is 
attached to an arbitrary raster choice which by nature, destroys the fundamental 2-dimensional 
nature of images. More sophisticated methods which aim at global optimization in some sense 
are very slow. 

SUMMARY OF THE INVENTION 

In view of the foregoing and other problems, disadvantages, and drawbacks of the 
conventional methods and structures, an object of the present invention is to provide a method 
and system which is much faster than the search of global optimum, but which still can deal with 
the 2-dimensional aspect of the images. 

In a first aspect, a method (and system) of digital printing includes forming a sequence of 
matrices P k with entries 0 or 1 where a 1 at some entry in some P k represents that this pixel will 
be printed at stage k, constructing as sequences of matrices I k with entries in [0, 1], so that I 0 = I, 
determining, when considering all pixels, a next pixel having a largest weight indicating that the 
next pixel is to be printed first, printing the pixel, and determining for each pixel of the 
remaining ones of pixels of the plurality of pixels a printing order of the remaining pixels such 
that subsequent pixels of the remaining ones of pixels of the plurality of pixels having a largest 
weight among the remaining pixels, are subsequently printed. 
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In a second Aspect, a method (and system) of printing includes forming a matrix of pixels, 
determining an order of printing of the pixels, the determining including finding a weight of the 
pixels and printing a pixel having a highest weight, and reordering the remaining pixels and 
printing a pixel having the greatest weight of the remaining pixels until all pixels have been 
printed. 

Further, a program storage medium is provided for storing program steps of the inventive 
methods. 

With the unique and unobvious features of the present invention, a method and system 
are provided which are much faster than the search of global optimum, but which still can deal 
with the 2-dimensional aspect of the images without treating each pixel independently of its 
neighbors. 



BRIEF DESCRIPTION OF THE DRAWINGS 



The foregoing and other objects, aspects and advantages will be better understood from 
the following detailed description of a preferred embodiment of the invention with reference to 
the drawings, in which: 

Figure 1 illustrates the use of dither masks for half-toning, according to conventional 
methods; 

Figure 2 shows a flow diagram which represents an exemplary method 200 of 

implementing the present invention; 
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Figure 3 shcAvs the structure of weights in a neighborhood V(i ? j) for a preferred 
embodiment of the present invention; 

Figure 4 illustrates an exemplary hardware/information handling system 400 for 
incorporating the present invention therein; and 

FIG. 5 illustrates a signal bearing medium 500 (i.e., storage medium) for storing steps of 
a program of a method according to the present invention. 

DETAILED DESCRIPTION OF A PREFERRED EMBODIMENT 
OF THE PRESENT INVENTION 

Turning to Figures 2-5, a detailed description of a preferred embodiment according to the 
present invention will be described. 

Prior to describing the preferred embodiment of the present invention and to make the 
description more precise, some notations will be introduced. Below, considered in detail is the 
case of digital printing with black and white dots of greyscale images. Other cases such as multi- 
tone and color will be discussed later. 

The original image will be denoted as the matrix I={I(i^ where 1 < i <Himage and 1 < 
j <Vimage. Each element I(i, j) of I is a real number between 0 and 1 , where 0 means "white" and 
1 means "black". 

The set of pixels to be printed after half-toning will be encoded in the matrix P = {P(i ; j)} 5 

where 1 < i < Himage and 1 < j <Vimage (i.e., this notation represents "horizontal" and 
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"vertical" coordinates of the image). Each element P(i ? j) of P is either 0 or 1, where again "0" 
means "white" and "1" means "black". It will be convenient to understand P as obtained from 
the constant 0 matrix by replacing as many 0s by 1 s as needed to represent I. 

Given an Himage-by-Vimage matrix M with entries in [0, 1], some neighborhood V(i„ j, 
M) is chosen for each (i, j), whose shape and size may either depend or not depend on (i 5 j) 5 and a 
set of weights associated to all pixels in V(i, j, M). Then, W(i, j, M) stands for the weighted 
average of the elements of M in V(i, j, M). 

Then, set GreyTotal= Zijl(ij). 

Hereinbelow and referring to Figure 2, a method (algorithm) 200 is described which 
constructs a sequence of matrices Pk with entries 0 or 1, the last one being the matrix P preferably 
sought by the present invention. In the process, there will be constructed sequences of matrices Ik 
with entries in [0,1], so that Io = I. 

INITIALIZATION: 

At step 201, first initialization is performed. That is, set k = 0 and define Po as the 
constant Himage-by-Vimage zero matrix. 

LOOP: 

In step 202, replace k by k+1 (i.e., as usual in languages such as FORTRAN, k now 
stands for what would have been k+1 before Step 1) and so forth for each iteration, and in step 

203 set Ik-Ik-i, P k =Pk-i. 
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In step 204, determine whether there has been a change in V(i, j, Mk). In steps 205 and 
206, compute all W(i, j, Ik) (weighted average) not computed so far and set W(i, j, Ik) = W(i, j, 

Ik-i) if W(i, j, Ik) is not computed again (i.e., when k = 1, one computes them all, and for k>l, one 
computes those which have been modified at the previous iteration). 
5 In step 207, order all W(i, j, Ik) by decreasing order, and in step 208 call Max(k) the 

bigger value of all W(i, j, Ik)s. 

In step 209, it is determined whether W(i, j, Ik) = Max(k). 

For all pairs (i, j) such that W(i, j, Ik) = Max(k) (i.e., a YES" in step 209), then in step 
^ 210, replace P k (i,j) = 0 with Pk(i,j)=l in P k , and I k (i,j) = I(i,j) withl k (i,j)-l(ij)-l inl k . Ifa 
i; jo "NO" results in step 209, then the process loops to step 21 1 in which nothing at (i, j) is changed. 

1 The process continues to step 212. In step 212, compute GreyTotal(k) = P k (i, j). 

j !: Then, if GreyTotal(k) < GreyTotal as determined in step 213, then the process loops to go to step 

202 for the next iteration, 
p Conversely, if GreyTotal < GreyTotal(k) as determined in step 214, and 

1 5 |GreyTotal-GreyTotal(k)| < |GreyTotal-GreyTotal(k-l)| as determined in step 215, then in step 

216setP-P k . 

Alternatively, if GreyTotal < GreyTotal(k) as determined in step 214, and 
|GreyTotal-GreyTotal(k)| > |GreyTotal-GreyTotal(k-l)| as determined in step 215, then in step 
217, set P = Pk-i. Thereafter, the process ends. 
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Of course, as would be known by one of ordinary skill in the art taking the present 
application as a whole, step 202 can easily be modified to vary the total number of black dots 
being printed, for instance to take account of the actual size of printed dots, compensate for 
originals which are too clear or too dark, just for a matter of taste, etc. 

To apply this method with a multi-tone printer, one prints the lightest grey darker than 
Max(k) (instead of necessarily a black dot in one bit per pixel printers), and replaces 

I k (iJ) = I(i,j)-l 

by 

Ik(i, j) = I(i, j) - "multibit grey level being printed" 
in steps 209-21 1 above. 

To apply this method to color images, the image I is decomposed as usual in the color 
planes accessible to the printer ( i.e., Cyan, Magenta, and Yellow, or Cyan, Magenta, Yellow, 
and Black) and then the method is used, as described above, in each plane. Then, in each color 
plane, either one bit or multi-tone printing can be accommodated readily. 

The invention can further be adapted to accommodate for the fact that some printers (e.g., 
such as some laser printers) require dark pixels to be grouped together, forming so-called half- 
tone dots, to compensate for isolated pixels that do not print in a predictable and trustable way. 

This accommodation is achieved by making the weights in V(i, j, M) depend also on what 
is printed at (i, j), for instance to favor printing near a pixel already designed to be printed. For 
instance, near a pixel already designed to be printed, the sum of the weights in V(i, j) may be 
greater than one, and otherwise may be smaller than one. 
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Thus; with the unique and unobvious features of the present invention, the tradeoff 
between performance and time required is optimized in that the time required by the inventive 
method, which may be greater than that for the conventional dither etc. techniques, is much 
smaller than for global optimization and other global methods. Conversely, the performance is 
5 much better than for the conventional dither etc. techniques. 

Figure 4 illustrates a typical hardware configuration of an information handling/computer 
system in accordance with the invention and which preferably has at least one processor or 
central processing unit (CPU) 411. 

The CPUs 41 1 are interconnected via a system bus 412 to a random access memory 
^10 (RAM) 414, read-only memory (ROM) 416, input/output (I/O) adapter 418 (for connecting 
^ peripheral devices such as disk units 421 and tape drives 440 to the bus 412), user interface 

"Ci adapter 422 (for connecting a keyboard 424, mouse 426, speaker 428, microphone 432, and/or 

s; other user interface device to the bus 412), a communication adapter 434 for connecting an 

IT information handling system to a data processing network, the Internet, an Intranet, a personal 

J:ijl5 area network (PAN), etc., and a display adapter 436 for connecting the bus 412 to a display 

device 438 and/or printer 439. As mentioned above, the printer 439 may be a digital printer or 
the like. 

In addition to the hardware/software environment described above, a different aspect of 
the invention includes a computer-implemented method for performing the above method. As an 
20 example, this method may be implemented in the particular environment discussed above. 
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1 r 

Such a method may be implemented, for example, by operating a computer, as embodied 
by a digital data processing apparatus, to execute a sequence of machine-readable instructions. 
These instructions may reside in various types of signal-bearing media. 

Thus, this aspect of the present invention is directed to a programmed product, 
comprising signal-bearing media tangibly embodying a program of machine-readable instructions 
executable by a digital data processor incorporating the CPU 41 1 and hardware above, to 
perform the method of the invention. 

This signal-bearing media may include, for example, a RAM contained within the CPU 
41 1, as represented by the fast-access storage for example. Alternatively, the instructions may be 
contained in another signal-bearing media, such as a magnetic data storage diskette 500 (Figure 
6), directly or indirectly accessible by the CPU 411. 

Whether contained in the diskette 500, the computer/CPU 41 1, or elsewhere, the 
instructions may be stored on a variety of machine-readable data storage media, such as DASD 
storage (e.g., a conventional "hard drive" or a RAID array), magnetic tape, electronic read-only 
memory (e.g., ROM, EPROM, or EEPROM), an optical storage device (e.g. CD-ROM, WORM, 
DVD, digital optical tape, etc.), paper "punch" cards, or other suitable signal-bearing media 
including transmission media such as digital and analog and communication links and wireless. 
In an illustrative embodiment of the invention, the machine-readable instructions may comprise 
software object code, compiled from a language such as "C", etc. 

With the unique and unobvious aspects of the present invention, a method and system are 

provided in which the tradeoff between performance and time required is optimized in that the 

time required by the inventive method, which is greater than that for the conventional error 
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diffusion, dithering etc. techniques, is much smaller than for global optimization and other global 
methods. Conversely, the performance is much better than for the conventional error diffusion, 
dithering etc. techniques. Indeed, the invention exploits between the two-dimensional aspect of 
the images than error diffusion methods and takes much more account to the image that dithering 
methods. Thus, the inventive method and system are much faster than the search of global 
optimum, but can still cope with the 2-dimensional aspect of the images. 

While a preferred embodiment of the present invention has been described above, it 
should be understood that it has been provided as an example only. Thus, those skilled in the art 
will recognize that the invention can be practiced with modification within the spirit and scope of 
the appended claims. 
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CLAIMS 

Having thus described our invention, what we claim as new and desire to secure by 
Letters Patent is as follows: 

1 . A method for digital printing represented as a matrix I, comprising: 

forming a sequence of matrices P k with entries 0 or 1 where a 1 at some entry in some P k 
represents that this pixel will be printed at stage k; 

constructing as sequences of matrices I k with entries in [0, 1], so that I 0 = I; 

determining, when considering all pixels in I k for all succesive values of k, a next pixel 
having a largest weight indicating that said next pixel is to be printed next , so that P k+i differs 
from P k by a zero at said next pixel in P k being replaced by a 1 at the same position in P k+1 ; 

printing said pixel; 

replacing the value of said pixel in I k by a 0 thus forming I k+1 ; and 

incrementing the value of k, until enough pixels have been printed to represent the overall 
darkness of I by the printed image. 

2. The method according to claim 1, wherein with each printing of a pixel, an order of printing of 
remaining pixels is redefined. 
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3 . The method according to claim 1 , further comprising reordering the pixels to be printed with 
each printing iteration. 

4. The method according to claim 1, wherein an original image is denoted as said matrix I = {I(i, 
j)}, where 1 < i <Himage and 1 < j <Vimage, and 

wherein each element I(i, j) of I is a real number between 0 and 1 where 0 represents 
"white", 1 represents "black" and intermediate values represent levels of grey. 

5. The method according to claim 4, wherein said matrix forming includes: 

setting k = 0 and defining P 0 as the constant Himage by Vimage zero matrix. 

6. The method according to claim 5, further comprising: 

with each iteration of redetermining, replacing k by k+1; and 
setting I k = I k _! and Pk = P k -i except for a single pixel. 

7. The method according to claim 6, further comprising: 

determining whether there has been a change in V(i, j, M k ); 

computing all weighted averages W(i, j, I k ) not computed in earlier iterations and setting 
W(i, j, Ik) = W(i, j, I k .i) if W(i, j, I k ) is not computed again, such that when k=l, all said weighted 
averages are computed, and for k>l, only said weighted averages which have been modified at a 
previous iteration are computed. 
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8. The method according to claim 7, further comprising: 

ordering all W(i, j, Ik) by decreasing order; and 
considering Max(k) as a larger value of all W(i, j ? Ik)s. 

9. The method according to claim 8 ? further comprising: 

determining whether W(i, j, Ik) = Max(k). 

10. The method according to claim 9, further comprising: 

for all pairs (i, j) such that W(i, j, Ik) = Max(k), replacing P k (i ? j) = 0 with P k (i, j) = 1 in 
P k , and I k (i 5 j) = I(i, j) with I k (i, j) = I(i, j)- 1 in I k . 

1 1 . The method according to claim 9 ? further comprising: 

computing GreyTotal(k) = Sj 0 - Pk(i, j). 

12. The method according to claim 1 1 ? wherein if GreyTotal(k) < GreyTotal, then a next 
iteration is begun, 

1 3 . The method according to claim 1 1 , further comprising: 

if GreyTotal < GreyTotal(k), and |GreyTotal-GreyTotal(k)| < |GreyTotal-GreyTotal(k-l)|, 

then set P=P k . 
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14. The method according to claim 1 1 , further comprising: 

if GreyTotal < GreyTotal(k), and |GreyTotal-GreyTotal(k)| > |GreyTotal-GreyTotal(k-l)|, 
then set P = P k .i. 

15. The method according to claim 1, wherein a total number of black dots being printed is 
variable. 

16. The method according to claim 15, further comprising: 

providing a multi-tone printer, such that a lightest grey darker than Max(k) is printed, and 
I k (i, j) = I(i, j) -1 is replaced by I k (i, j) = I(i, j) which represents a multibit grey level being 
printed. 

17. The method according to claim 9, further comprising: 

compensating for isolated pixels by making weights in V(i, j, M) depend on what is 
printed at (i, j). 

"18. The method according to claim 1 , further comprising: 

encoding the set of pixels to be printed after half-toning, in the matrix P={P(i, j)}, where 
1 < i < Himage and 1 < j <Vimage, wherein each element P(ij) of P has a value of either 0 or 1, 
where "0" represents "white" and "1" represents "black". 
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19. The method according to claim 1 8, further comprising: 

given a Himage by Vimage matrix M with entries in [0,1], selecting a neighborhood V(i, 
j, M) for each (i, j), whose shape and size selectively depends or not on (i, j), and a set of weights 
associated to all pixels in V(i, j, M). 

20. The method according to claim 19, wherein W(i, j, M) represents a weighted average of the 
elements of M in V(i, j, M), and GreyTotal= Xjj I(i, j). 

21. A method of printing, comprising: 

forming a matrix of pixels; 

determining an order of printing of said pixels, said determining including finding a 
weight of said pixels and printing a pixel having a highest weight; and 

reordering the remaining pixels and printing a pixel having the greatest weight of the 
remaining pixels until all pixels have been printed. 

22. The method according to claim 21, wherein the pixels are printed in turn based on the 
darkness of the local image being printed. 

23. A method for digital printing, comprising: 

forming a sequence of matrices Pk with entries 0 or 1 where a 1 at some entry in some Pk 
represents that this pixel will be printed at stage k; 
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constructing as sequences of matrices I k with entries in [0, 1], so that Io = I; 
determining, for a plurality of pixels in I k for all successive values of k, a next pixel 
having a largest weight indicating that said first pixel is to be printed next; 
printing said pixel; and 

determining for each pixel of the remaining ones of pixels of said plurality of pixels a 
printing order of said remaining pixels such that subsequent pixels of said remaining ones of 
pixels of said plurality of pixels having a largest weight among the remaining pixels, are 
subsequently printed. 

24. A system for printing, comprising: 

means for forming a matrix of pixels; 

means for determining an order of printing of said pixels, said determining including 
finding a weight of said pixels and printing a pixel having a highest weight; and 

means for reordering the remaining pixels and printing a pixel having the greatest weight 
of the remaining pixels until all pixels have been printed. 

25. The system according to claim 24, wherein with each printing of a pixel, an order of printing 
of remaining pixels is redefined. 

26. The system according to claim 24, wherein said reordering means reorders the pixels to be 
printed with each printing iteration. 
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27. The system according to claim 24, wherein an original image is denoted as said matrix I = 
{I(i, where 1 < i <Himage and 1 ^ j <Vimage, and 

wherein each element I(i, j) of I is a real number between 0 and 1 where 0 represents 
"white", 1 represents "black" and intermediate values represent levels of grey. 

28. The system according to claim 27, wherein said matrix forming means includes: 

means for setting k = 0 and defining Po as the constant Himage by Vimage zero matrix. 

29. The system according to claim 28, further comprising: 

with each iteration of redetermining, means for replacing k by k+1 ; and 
means for setting Ik = Ik-i and Pk= Pk-i except for a single pixel. 

30. The system according to claim 29, further comprising: 

means for determining whether there has been a change in V(i, j, Mk); and 
means for computing all weighted averages W(i, j, Ik) not computed in earlier iterations 
and setting W(i, j, Ik) = W(i, j, Ik-i) if W(i, j, Ik) is not computed again, such that when k=l , all 
said weighted averages are computed, and for k>l, only said weighted averages which have been 
modified at a previous iteration are computed. 

31. The system according to claim 30, further comprising: 
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means for ordering all W(i, j, I k ) by decreasing order; and 
means for considering Max(k) as a larger value of all W(i, j, I k )s. 

32. The system according to claim 31, further comprising: 

means for detennining whether W(i, j, I k ) = Max(k). 

33. The system according to claim 32, further comprising: 

for all pairs (i, j) such that W(i, j, I k ) = Max(k), means for replacing P k (i, j) = 0 with P k (i, 

j) = 1 in P k , and I k (i, j) = I(i, j) with I k (i, j) = I(i, j)-l in Ik- 

34. The system according to claim 32, further comprising: 

means for computing GreyTotal(k) = Ey P k (i, j). 

35. The system according to claim 34, wherein if GreyTotal(k) < GreyTotal, then a next iteration 
is begun. 

36. The system according to claim 34, further comprising: 

if GreyTotal < GreyTotal(k), and |GreyTotal-GreyTotal(k)| < |GreyTotal-GreyTotal(k-l)|, 

means for setting P=P k . 
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37. The system according to claim 34, further comprising: 

if GreyTotal < GreyTotal(k), and |GreyTotal-GreyTotal(k)| > |GreyTotal-GreyTotal(k-l)|, 

means for setting P = Pk-i- 

38. The system according to claim 24, wherein a total number of black dots being printed is 
variable. 

39. The system according to claim 38, further comprising: 

a multi-tone printer for printing, such that a lightest grey darker than Max(k) is printed, 
and I k (i, j) = I(i, j) -1 is replaced by I k (i, j) = I(i, j) which represents a multibit grey level being 
printed. 

40. The system according to claim 32, further comprising: 

means for compensating for isolated pixels by making the weights in V(i, j, M) depend on 

what is printed at (i, j). 

41. The system according to claim 24, further comprising: 

means for encoding the set of pixels to be printed after half-toning, in the matrix P={P(i, 
j)}, where 1 < i < Himage and 1 < j <Vimage, wherein each element P(ij) of P has a value of 
either 0 or 1, where "0" represents "white" and "1" represents "black". 

42. The system according to claim 4 1 , further comprising: 
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given a Hirriage by Vimage matrix M with entries in [0,1 ], means for selecting a 
neighborhood V(i, j, M) for each (i, j), whose shape and size selectively depends or not on (i, j), 
and a set of weights associated to all pixels in V(i, j, M). 

43 . The system according to claim 42, wherein W(i, j , M) represents a weighted average of the 
elements of M in V(i, j, M), and GreyTotal= Sy I(i, j). 

44. The system according to claim 24, wherein the pixels are printed in turn based on the 
darkness of the local image being printed. 

45 . A signal-bearing medium tangibly embodying a program of machine-readable 
instructions executable by a digital processing apparatus to perform a method of printing, said 
method comprising: 

forming a matrix of pixels; 

determining an order of printing of said pixels, said determining including finding a 
weight of said pixels and printing a pixel having a highest weight; and 

reordering the remaining pixels and printing a pixel having the greatest weight of the 
remaining pixels until all pixels have been printed. 
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46. A signal-bearing medium tangibly embodying a program of machine-readable instructions 
executable by a digital processing apparatus to perform a method of printing, said method 
comprising: 

forming a sequence of matrices Pk with entries 0 or 1 where a 1 at some entry in some Pk 
means that this pixel will be printed at stage k; 

constructing as sequences of matrices Ik with entries in [0 ? 1], so that Io= I; 

determining, when considering all pixels in Ik for all successive values of k, a next pixel 
having a largest weight indicating that said next pixel is to be printed next , so that Pk+i differs 
from Pk by a zero at said next pixel in Pk being replaced by a 1 at the same position in Pk+i; 

printing said pixel; 

replacing the value of said pixel in Ik by a 0 thus forming Ik+i; and 

incrementing the value of k, until enough pixels have been printed to represent the overall 
darkness of I by the printed image. 
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METHOD AND SYSTEM FOR DIGITAL PRINTING BY DENSITY SCHEDULING 

ABSTRACT OF THE DISCLOSURE 

A method (and system) of printing, includes forming a matrix of pixels, determining an 
order of printing of the pixels, the determining including finding a weight of the pixels and 
printing a pixel having a highest weight, and reordering the remaining pixels and printing a pixel 
having the greatest weight of the remaining pixels until all pixels have been printed. 
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